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Infra-red stability of Yukawa and soft-breaking fixed points 
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We investigate the infra-red stability of the fixed points in the evolution of the Yukawa 
couplings, A-parameters and soft scalar masses in a broad class of supersymmetric theories. 
We show that the issue of stability is essentially determined in all three cases by the 
eigenvalues of the same matrix. In a very wide range of physically interesting theories it 
follows that, in the asymptotically free case, the existence of stable infra-red fixed points 
for the Yukawa couplings implies stable infra-red fixed points for the A-parameters and 
soft scalar masses. 
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The predictive power of the supersymmetric standard model and its extensions may 
be enhanced if the renormalisation group (RG) running of the parameters is dominated by 
infra-red (IR) stable fixed points. Typically these fixed points are for ratios; for instance, 
of the Yukawa coupling to the gauge coupling^, or the A-parameter to the gaugino mass. 
Moreover, even if these ratios have not attained their fixed point values at the weak scale 
(which in fact is usually the case), the couplings may be determined by quasi-fixed point 
behaviour @ . Here the value of the coupling at the weak scale is independent of its value at 
the unification scale. For this scenario, the Yukawa couplings at the unification scale need 
to be large. There is a considerable literature devoted to the consequences of fixed-point 
and quasi-fixed point behaviour in the standard model and the MSSM, and in extensions 
thereof p|— [53]. A necessary condition for the success of both the fixed point and the 
quasi-fixed point scenarios is the existence of stable infra-red fixed points. Hence it is of 
interest to be able to determine for a given theory the nature of the infra-red fixed points, 
and this is the question we address in this paper. We consider N = 1 supersymmetric 
gauge theories, and investigate the IR stability of the fixed point structure of the Yukawa 
couplings, the soft-breaking A-parameters and the soft-breaking masses^. The stability of 
the fixed points in each case is determined by the positivity of the eigenvalues of a matrix 
M—pQI where Q is the coefficient of the one-loop gauge /3-function, and the matrix M and 
the parameter p depend on the type of coupling. Our main result is that for a theory for 
which the anomalous dimension matrix for the matter fields is diagonal, there are matrices 
R, S (which are r x s and s x r respectively) and D (which is r x r and diagonal) such 
that the matrix M is given by DRS for the Yukawa case, RSD for the A-parameter case 
and SDR for the soft mass case. It is easy to show that the non-zero eigenvalues of the 
three matrices DRS, RSD and SDR coincide, and hence the eigenvalues determining the 
stability of the fixed points for all the couplings in the theory are given in terms of one 
basic set. Moreover in a more restricted class of theories (but which still includes all cases 
commonly considered) we can show that these non-zero eigenvalues are all positive. This 
implies that in the case of negative Q, all the fixed points are infra-red stable. 

For aniV = l supersymmetric gauge theory with gauge group H a G a , superpotential 

w = iy^^ + ^ y $<$ is (i) 



1 often called Pendleton-Ross (PR) fixed points |T[ 

2 Stability of the Yukawa couplings alone has been considered for a range of theories in Ref. |] 
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and soft breaking terms given by 



£sb = (m 2 )2>Vj + (^h^fofofa + + \Y,M a \ a \ a + h.c.j 



(2) 



(where is the scalar component of $, and A is the gaugino), the /^-functions for Y, h, b 
and m 2 are given by 



167r 2 4 1)<ifc 



Y^ p P k p + (k^i) + (k^j), 

h ijiph l + Y i i l X k l + {k^i) + {k^ j), 



16n 2 [3l mj) = b il P k l + ^Y^ l Y lmn b mn + fj, il X^, 



16tt 2 [/3, 



l -Y ipq Y^ n {m 2 y n + 1 -Y^Y pqn {m 2 r i + Tf ipq Y^{m 2 )\ 



+ h ipq V™ - SMaMZglCtfj + 2gi(R a A y j Tr[R a A m 



a _2n 



(3a) 
(36) 
(3c) 

(3d) 



where 



(4) 



pi . _ ly-ikly _ 9 2mri . 
3 2 Y 3kl Z 9a^i 3i 

X i j = W kl Y jkl + AMglC^S i j . 

Here g a is the gauge coupling for the gauge group G a and Cf is the quadratic Casimir for the 
representation R A of <frj. The matrix P is related to the one- loop anomalous dimension by 
' = j^P. Now suppose that the distinct, independent couplings in the superpotential 
are denoted Y a , a = 1, . . . r, and a = 1, . . .q, and that the distinct group multiplets 
are denoted I = l,...s. Let us assume that P l j is diagonal, with = -P/O^ for 
all $i in Further, let us assume that one of the couplings g a (where G a ^ U(l)) and 
its corresponding gaugino mass M a are dominant in the RG evolution. Henceforth we 
denote these parameters by g, M and suppress the remaining gauge couplings and gaugino 
masses. We may write 



IaUc 



2g 2 d 



(5) 



where y a = Y 2 . Here Si a is related to the dimensionality of the $7; clearly it is zero unless 
the superpotential term involving y a contains $7. The one- loop RG equation for Y a is 

d 



dt 



Y a — Ya ^ RalPl 



(6) 



where 16n 2 t = In/i and R a i takes the value 0, 1, 2 or 3 according as the superpotential 
term involving y a contains 0, 1, 2 or 3 $/s. Note that the zero entries of R a i coincide 
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with those of Si a , and we have R a j = 3 for all a. As an example, for the usual MSSM 
superpotential, retaining only 3rd generation Yukawas and (73, 

W = Y t QtH 2 + Y h QbH x + Y t LtH u (7) 

we have 



P Q = Y t 2 + Y b 2 -lgl P L = Y T 2 , P t = 2Y t 2 -lgi 
P b = 2Y 2 -\gl P T =2Y 2 1 Ph 2 =ZY?, P Hl =?,Y 2 + Yl 



(8) 



so that R Yt Q = Ry b Q = 1, Ry t t = Rv b b = 1, and S QYt = S QYb = 1, S tYt = S bYb = 2 etc. 

For the moment we shall focus on the parameters Y^ k , h 1 ^ and m 2 ; we shall discuss 
the \i and b parameters later. We assume that the distinct independent trilinear soft- 
breaking couplings are in one-to-one correspondence with the Y a , and we define A a = 
We then find 

16tt 2 ^ = 2 R a iS ip ypAp + 4 ^ g 2 MR ctI C I . (9) 



1,(3 I 

Similarly, we find, assuming a common soft mass mj for each multiplet 

a, J a 

The one-loop RG equations for g and M are 



16tt 2 ^ = 2 S Ia y a Rajm 2 j + 2 ^ S Ict y a A a A* a - 8g 2 MM*d. (10) 



| 5 = Qg 3 , jM = 2Qg 2 M, (11) 
and it follows from Eqs. @, (pTTf ) that the RG evolution of y a = %, A a = ^ and 



^y a = 2y a ^^i? Q /P/ - g 2 Q 



d A a = 2g 2 R a iS ip ypAp - 2Qg 2 A a + 4 ^ g 2 R aI C I (12) 



dt 

i,p 1 
<] ~ j = 2g 2 J2 Si a y a R a jrh 2 - 4Qg 2 mj + 2g 2 £ S Ia y a A a A* a - 8g 2 C T . 



—m 

dt 

a, J 



These equations may be rewritten in the more compact form 

g- 2 j t y = 2D{RSy-U) 



g- 2 ^-A = 2(RSD -QI)A + 4RC (13) 

g- 2 ^-fh 2 = 2(SDR - 2QI)m 2 + 2SDC - 8C, 
dt 



where R, S and D are r x s, s x r and r x r matrices with D a p = y a 5 a p, and the column 
vectors U and C are defined by 

E/* = 2(RC) a + Q, C a =A a A* a . (14) 

We now see the fixed point structure quite clearly. Turning firstly to the Yukawa couplings, 
y a = is always a possible fixed point for each a. The most general fixed point satisfies 

y a = (a e J) 

(15) 

(RSy-U) a =0, (ae{l,...r}\J) 

where J is an arbitrary subset of {1, . . . r}. If this condition is satisfied then we also see 
that 

A a = -1 (a e {1, . . .r}\J) 



A 



a 



^ ^2(RC) a - J2( RSD )^ (« e J ) 



(16) 



represents a fixed point. (If Q = and J ^0, then there are in general no fixed points 
for A a with a G J.) The existence of this fixed point, corresponding to a universal 
A-parameter, A a = — M (in the case J =0), was first remarked in Ref. (and the 
importance of this value for A in the special case of a finite theory was realised earlier, in 
Ref. ||25|| ). If (SDR — 2QI) is invertible then the fixed point for rh 2 is given by 

m 2 = (SDR- 2QI)' 1 {AC - SDC), (17) 



while otherwise there is no obvious closed form expression. It is easy to show from Eq. fll3| ) 
that in the case J =0, the quantity (Rrh 2 ) a has the fixed point (Rrh 2 ) a = 1 for all a. This 
"sum rule" has been investigated in Refs. . We should mention also that in the special 
case where the Yukawa fixed points satisfy (Sy — 2C)i = for all / (corresponding to 



the so-called P = case[^7|) then it is straightforward to show that m 2 = | is a fixed 



point for all /. 

In any event, our main concern here is with the IR stability of the fixed points rather 
than their actual values. The stability of the fixed points of the evolution equations for 
a generic set of parameters Aj is determined by linearising the evolution equations about 
the fixed points A* so that they adopt the form -^SXi = ffrl^.^Aj, where 5\ = A — A*. 
The system is IR stable if all the eigenvalues of ff 1 1 A » have positive real parts. In the case 
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at hand the task is fairly simple, the soft-breaking equations being linear already, and the 
only subtlety arises for the Yukawa coupling case. We find 

^ = 28 af3 (RSy - U) p + (DRS) aP . (18) 

Let us first of all consider the simplest case, of a fixed point with all of the y a non-zero, in 
other words J =0. The Yukawa fixed points then simply satisfy RSy — U = and it follows 
from Eq. ([TJ) that the Yukawa stability matrix is just DRS. We then see immediately 
from Eq. (|T3| ) that the stability matrices for the A and m 2 are RSD — QI and SDR — 2QI. 
Hence for IR stability we require DRS, RSD — QI and SDR — 2QI to have positive 
eigenvalues. The appearance of combinations of R, S and D in each of these stability 
matrices seems quite remarkable. Even more interestingly, the three matrices DRS, RSD 
and SDR have the same set of non-zero eigenvalues. This follows from the easily-proven 
result that MN and NM have the same non-zero eigenvalues for any m x n and n x m 
matrices M, N. 

We have some immediate and surprisingly powerful results. In an asymptotically free 
theory (Q < 0) such that the Yukawa couplings flow to non-zero PR fixed points, the 
A-parameters and soft masses also flow to fixed points. A non-asymptotically-free theory 
(Q > 0) cannot have completely stable fixed points for both A-parameters and soft masses 
unless the number of Yukawa couplings and the number of fields is the same. This is 
because unless the matrices R and S are square, then either RSD or SDR will have some 
zero eigenvalues (except possibly in pathological cases where, RSD (say) has fewer than r 
independent eigenvectors). 

Now let us turn to the case of a fixed point where some of the y a (p, say) are chosen 
to be zero. Let us re-order the y a so that J = {r — p + 1, . . . r}, and let us partition the 
y into y a , a = 1, . . . r — p and a = r — p + 1, . . . r. It follows from Eq. ( |T5| ) that the 
stability matrix for the Yukawa couplings now takes the block form 



T 3 



(19) 



where 

T laf3 = (DRS) af3 , T 2a - p = {DRS) a - p , T^p = ^{RSy - U)^ (20) 

and (with some abuse of notation) R, S and D represent the relevant matrices after re- 
ordering the couplings. Hence for stability of the Yukawa fixed points we require {RSy — 
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U) q > for [3 = r — p + 1, . . .r together with the positivity of the eigenvalues of the matrix 
T\. The stability matrices for the A and m 2 are still given by RSD — QI and SDR — 2QI, 
and it still follows that the non-zero eigenvalues of DRS, RSD and SDR coincide-though 
there are now additional zero eigenvalues, since D contains p zeroes on the diagonal. The 
eigenvalues of RSD — QI consist of — Q (p times) together with the eigenvalues of T\ — QI. 
Hence this fixed point can only be stable in the asymptotically free case. 

Finally we should say a word about the /i and b parameters. The evolution of the \i 
parameter is governed by 

9~ 2 f f ^ = V&(RSy - 2RC) & , (21) 

where R^i takes the value 0, 1 or 2 according as the \x term contains 0, 1 or 2 $s. The 
only fixed point for [i (or any combination such as ^, -j^) is \ia = for each a (except in 
the special case of a one- loop finite theory with P = 0, for which any value of \i is a fixed 
point). The stability of = requires (RSy — 2RC)et > for each a. In the case of a 
theory with no gauge singlets, we have Yi mn b rnn = in Eq. (0). In that case we have 



g - 2 ^-B = -2QB + 2RSDA + ARC, (22) 

where = ^ & M - Eq. (^2|) clearly determines the fixed point for B in terms of that for 
A. Moreover, the stability of the fixed point for B^ simply requires that Q be negative. 
However, in a theory with gauge singlets, the stability of B& is determined by a new matrix 
which has no obvious connection with those which have arisen hitherto. 

Let us illustrate our results with some simple examples. We begin with a simplified 
version of the NMSSM, with the superpotentiali 

W = Y t QtH 2 + \NH±H 2 - fiV 3 . (23) 

The one-loop anomalous dimensions are given by 



P Q = Y t 2 -lgl P t = 2Y t 2 -lgl P H2 =3Y t 2 + X' 



? Hl = X\ P N = 2(A 2 + k 2 ' 



(24) 



3 Here (and also in the next example of the MSSM) we ignore all except the third generation 
Yt coupling. It is straightforward to generalise to the case of diagonal couplings for the other 
generations. 
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Writing Y t = etc, we find that our matrices R : <§, D are given in the NMSSM case by 



1110 
R= | 1 1 1 

3 



S 



(1 0^ 

2 

3 1 
1 

\0 2 2/ 




(25) 



The Yukawa coupling evolution in this model was analysed in Ref. | JL8|| . The only non- 
trivial fixed point with Y 2 , A 2 , k 2 > is 

A = fc = 0, Y 2 = ^ 



and is clearly stable, since the stability matrix of Eq. ([19]) is just 
Eq. ([1^) we find the fixed point for the A to be 




A 



Y t 



•1, A\ — jg, A 



0, 



and from Eq. (|T^) we find the fixed points for the m 2 to be 



From 



(26) 



41 
54' 



11 

27' 



H 2 



18' 



(27) 



in agreement with Ref. [JT9J. Again these fixed points for the A and m 2 are clearly stable, 
since Q is negative and T\ = |. Our formalism excludes the possibility of determining 
fixed point values for A-parameters corresponding to mixing between the generations. 

For our second example we take the model with the MSSM superpotential of Eq. (|7|). 
The matrices R and S are given (ordering the fields as Q, L, t, 6, r, i?2, Hi) by 



10 10 10 

S =(1001001 

10 10 1 



s 



f 1 


1 










1 


2 











2 











2 


3 








Vo 


3 


1/ 



y t 2 

Y 2 






Y 2 



(28) 



There is no fixed point with Y 2 > 0. However, in the case where we take Y T = and Y^, 
Y"b non-zero at the fixed point, we find a fixed point with 



We find in Eq. (|19|) that at the fixed point T\ 



2 ±\ 

i ^ J and T3 = 1. The non- vanishing 

eigenvalues of T\ are | and |. Since moreover Q < for the MSSM, we deduce that the 
fixed points are stable. Evidently the above fixed point values for Y tl Yi correspond to 
large tan/3 (specifically tan/3 = mt/rrib); the RG evolution of the Yukawa couplings in this 
region was studied in Ref. |10|1 . It would be interesting to extend this analysis to include 
the soft-breaking terms. 

For our third example we take a model with superpotential 



W 



YiiQiti~H.ii. 



(30) 



We have 



Rn,i 

SiAI 



1, 
1, 



RiI,I 

Siai 
Duai 



R, 



1, 

2, Suai 



IAI 



= 1, 

3, 



(31) 



yu, 



with all other entries in R, S and D being zero. The Yukawa fixed point is yn = y, where 
V = (^rqrfyj f° r an I (assuming we take all the y non-zero). In this case the fixed point 
for every A is An = — 1, and the fixed points for m 2 are 



^2 



il 



3y 
Q ' 



(32) 



For stability it is sufficient to consider RS. It is straightforward to compute the eigenvalues 
explicitly in this case. We find that RS — XI is the n\ 



x n g a matrix given by 



([h-\]I + K 21 21 

21 [5-\}I + K 21 



RS-XI 



V 



21 



21 



21 



21 
21 



21 
21 



(33) 



21 [5-\}I + K/ 



where K is an n g x n g matrix with every entry 1, and / is the n g x n g unit matrix. By 
subtracting the last column from every other column, and then adding every row to the 
last row, we conclude that 



det(RS - XI) = det[(3 - X)I + K] nB ~ x det[(3 + 2n g - X)I + K]. 
By similar operations we can show that 

det[xl + K] = x n9 ~ 1 (x + n g ), 



(34) 



(35) 
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and it follows that the eigenvalues of RS are 3 ([n g — l] 2 times), n g + 3 ([n g — 1] times), 
2n g + 3 ([rig — 1] times) and 3n g + 3 (once). These are all positive. Hence again we 
find stability in the asymptotically free case. It is easy to combine our second and third 
examples to obtain the n g generation theory which was considered as a model for strong 
unification in Ref. [ESI. 

As our final example we will take a semi-realistic GUT with gauge group SU3 ® SU3 (g> 
SU3, and a matter content consisting of n sets each of the representations X = (3,3, 1), 
Y = (1, 3, 3) and Z = (3, 1, 3). The superpotential for the theory is : 



W 



1 
3! 



(A1X 3 + X 2 Y 3 + A 3 Z 3 ) + pXYZ. 



(36) 



Here A1X 3 = (Ai) ^ X a X ( gX 7 , where a, /3 • ■ ■ = 1 • ■ ■ n. If we set the three gauge couplings 
all equal to g then it is easy to see that they remain equal under renormalisation. In what 
follows we will suppose that we have (A 2 )^ = (A 2 )5^, and (p 2 )p = p 2 8p. In this case we 
have (note that this is a case with s > r) 



R 



f 3 




Vi 








( 2 












3 

1 


s= 


2 







D = 


3 

J 


\o 





2 




V 










A! 








\2 

A 3 










(37) 



At the fixed point we have 6A 2 + 9p 2 = Q + 16, and hence we must take Q > —16 to 
ensure physical fixed point values. We find the non-zero eigenvalues of DRS etc to be 
Q + 16 — 9p 2 (twice) and Q + 16, which are positive provided the A 2 and p 2 are positive. 



The implications for the stability of the infra-red fixed points were discussed in Ref. |12 
Again stability is assured in the case Q < 0. 

Each example we investigated had the property that the non-zero eigenvalues of DRS 
were positive (for physical fixed-point values for the y, in other words each diagonal entry 
of D non-negative at the fixed point) leading to IR stability of the fixed points for A, m 2 
for an asymptotically free theory. The obvious question is whether this is true in general. 
We have so far been unable to answer this question completely; however, as a first step we 
can show that if S may be written S = D\R T D 2l where D\ and D 2 are diagonal matrices 
with positive entries along the diagonal, and if D has non-negative entries, then DRS has 
non-negative eigenvalues. All the examples we have considered have this property. For 
instance, in our first example, Di = diag(l, 2, 3, 3, 6), D 2 = diag(l, in our second 

example, Di = diag(l, 3, 2, 2, 6, 3, 3), D 2 = diag(l, 1, -|); in our third example, Du^ = 1, 
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Duj = 2, Duj^ij = 3, D 2 = /; and in our fourth example, Di = |/, D 2 = diag(l, 1, |). 
It is possible to write down matrices R and S, subject to the constraints that R a i = 3 
and that R corresponds to a theory with P diagonal, which do not have the property 
above; but they do not correspond to any physically interesting model, and in any case we 
have not been able to construct any examples which lead to non-zero eigenvalues of DRS 
with negative real parts. It is tempting to speculate that this is a general result. 

We thus conclude that the focusing of Yukawa couplings and soft breakings observed 
for example, in Ref. |19|] for the MSSM is not in fact specific to that theory but a general 
phenomenon; the existence of stable infra-red fixed points for the Yukawa couplings implies 
(given asymptotic freedom) stable infra-red fixed points for the A-parameters and soft 
scalar masses. In Ref. [[12) we proposed that universality of soft masses and couplings at 
the unification scale might be associated with IR fixed points. Our demonstration here 
that there exists a strong relationship between the IR stability of the fixed points for these 
masses and couplings with the stability of the fixed points in the Yukawa sector makes this 
scenario even more plausible. 
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